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Abstract 

An n-simplex is said to be n-well-centered if its circumcenter lies in its interior. We introduce 
several other geometric conditions and an algebraic condition that can be used to determine 
whether a simplex is n-well-centered. These conditions, together with some other observations, 
arc used to describe restrictions on the local combinatorial structure of simplicial meshes in 
which every simplex is well-centered. In particular, it is shown that in a 3-well-centered (2- 
well-ccntcrcd) tctrahedral mesh there are at least 7 (9) edges incident to each interior vertex, 
and these bounds are sharp. Moreover, it is shown that, in stark contrast to the 2-dimensional 
analog, where there are exactly two vertex links that prevent a wcll-ccntcrcd triangle mesh in R 2 , 
there are infinitely many vertex links that prohibit a well-centered tetrahedral mesh in R 3 . 

1 Introduction 

An n-dimensional simplex is n-well-centered if its circumcenter lies in its interior. More gen- 
erally, it is k -well- centered if each of its fc-dimensional faces is k- well-centered. It is completely 
well-centered if it is fc-well-centered for each k, 1 < k < n [15]. Several authors have noted 
the possible application of well-centered meshes to particular problems [15]. Among these are 
Nicolaides [10] and Sazonov et al. [11, 12], who have discussed the covolume method and its 
application in electromagnetics simulations. Also, Kimmcl and Scthian [5] described an algo- 
rithm for numerically solving the Eikonal equation on triangulated domains. Their algorithm, 
which can be used to compute geodesic paths on triangulated surfaces, is described first for 
acute triangulations (i.e., 2- well-centered triangulations) and requires additional work for tri- 
angulations that are nonacute. Ungor and Sheffer [13] used acute planar triangulations when 
they introduced the tent-pitching algorithm for space-time meshing. Well-centered meshes also 
find application within Discrete Exterior Calculus (DEC), a framework for designing numerical 
methods for partial differential equations [4, 3]. In DEC, a sufficient (but not necessary) condi- 
tion for discretizing the Hodge star operator so that it is represented by a diagonal matrix is to 
use a well-centered mesh. The diagonal matrix leads to efficient numerical solution. 

Constructing well-centered meshes is a nontrivial task, so researchers have put effort into 
finding ways to work around an algorithmic requirement for well-centered meshes. Recently, 
however, there has been progress towards making well-centered meshes through mesh optimiza- 
tion. In [16] we described a heuristic for obtaining well-centered meshes of planar domains by 
starting with a given triangle mesh of the domain and relocating the interior vertices of the mesh 
to minimize a cost function defined on the coordinates of those vertices. We later generalized 
the cost function, making it possible to optimize meshes of any dimension [17]. In [15] we give 
a variety of examples of tetrahedra that are and are not well-centered, and we show that it is 
possible to mesh many simple shapes in R 3 with completely well-centered tetrahedra. In many 
cases, the completely well-centered tctrahedral meshes of [15] were obtained by creating initial 
relatively high-quality meshes by hand and applying the optimization method to the initial 
meshes. 



Not all simplicial meshes can be made well-centered by optimizing the cost functions of [16] 
and [17]. In some cases the combinatorial properties of the mesh prevent the mesh from be- 
coming well-centered. In such cases the mesh will not be well-centered for any choice of vertex 
coordinates. This paper develops theory and intuition related to the geometry of well-centered 
simplices and applies the geometric conditions to investigate combinatorial properties of well- 
centered meshes. 

Before discussing the specific results of this paper, a few comments about terminology are 
in order. The term well-centered may be used without a qualifying dimension to refer to the 
general concept, or may refer to one of the more precise terms if the context makes clear which 
more precise term is appropriate. When speaking of a triangle, for example, well-centered simul- 
taneously means 2-well-centered and completely well-centered, since every simplex is trivially 
1-well-ccntered. All of these definitions can be applied to simplicial meshes or simplicial com- 
plexes embedded in Euclidean space. Thus if a simplicial mesh is said to be fc-well-centered, 
this means that every /c-dimensional simplex appearing in the mesh properly contains its cir- 
cumcenter. 

2 Results 

After giving definitions and notation in Sec. 3, we introduce in Sec. 4 several geometric and 
algebraic conditions for an n-simplex to be well-centered. Sections 5 and 6 apply the theory 
developed in Sec. 4 to establish conditions on the combinatorial structure of the neighborhood 
of a vertex in a well-centered tetrahedral mesh. Finally, Sec. 7 records some observations specific 
to constructing well-centered meshes of the cube, and Sec. 8 offers some concluding thoughts. 
We enumerate the contributions of this paper in more detail in the following paragraphs. 

In Sec. 4 we prove three new conditions on when an n-simplex is n-wcll-centered, each phrased 
in terms of the location of a vertex Vi given the facet of the simplex opposite Uj. These conditions 
are (a) a necessary condition expressed in terms of geometry (Prop. 3 — the Cylinder Condition), 
(b) a sufficient condition expressed geometrically (Prop. 8 — the Prism Condition) , and (c) a both 
necessary and sufficient condition expressed in terms of cubic polynomial inequalities (Prop. 11). 
The two geometric conditions are generalizations to higher dimensions of conditions in R 2 . 

Section 5 investigates combinatorial properties that follow from the results in Sec. 4. (d) We 
prove a new combinatorial condition that must be satisfied by the link of an interior vertex 
in an n-well-ccntcrcd mesh in R™ (Thm. 13). (c) As an easy corollary we show that in a 3- 
well-centered tetrahedral mesh in R 3 , every interior vertex has at least seven incident edges 
(Cor. 14). (f) We show that, in stark contrast to the analogous case in R 2 , where there are only 
two vertex links that cannot appear in a 2-well-centered mesh, there are infinitely many vertex 
links that cannot appear in a 3-well-centered mesh in R 3 (Cor. 15). (g) We also construct an 
infinite family of vertex links that can appear in a completely well-centered tetrahedral mesh 
in R 3 . (h) The section closes by showing that if a vertex link can appear in a 3-well-centered 
mesh and the vertex link satisfies some minor additional conditions, then degree three vertices 
can be successively inserted into the vertex link to create an infinite family of vertex links that 
can appear in a 3-well-centered mesh (Prop. 18). 

Section 6 develops combinatorial conditions that 2-well-centered tetrahedral meshes in R 3 
must satisfy, (i) We prove in Thm. 22 that no triangulation of S 2 on m vertices with a vertex 
of degree at least m — 3 can appear in a 2-well-centered tetrahedral mesh in R 3 . (j) It follows 
that in a 2-well-centered tetrahedral mesh in R 3 , every interior vertex has at least nine incident 
edges (Cor. 23). (k) We show that vertices of degree three can be inserted into or deleted from 
triangulations that permit 2-well-centercd neighborhoods to create other triangulations that 
permit 2-wcll-ccntcrcd neighborhoods (Prop. 24). (1) Vertices of degree four can also be added 
to such triangulations (Prop. 25). 

At several points in the paper, we make claims about additional results beyond what is 
actually proved in this paper. Further details about some of these claims appear in the first 
author's dissertation [14]. 
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3 Definitions and Notation 

We begin by introducing some definitions and notation that will be used throughout the paper. 
A simplex is referred to with a Greek letter, usually a or r. A superscript for a simplex indicates 
the dimension, so, for example, a n is an n-simplcx. The notation a n = [vqVi . . . v n ] is used to 
indicates that a n is the convex hull of the n + 1 vertices vo, v\ , . . . , v n . It is assumed that the 
vertices of a simplex are in general position, i.e., that the vertices are affinely independent, so 
a n is fully n-dimcnsional. The circumcenter of a simplex a is denoted c(tr). For an n-simplcx 
a n embedded in R n , c(<r n ) is the unique point which has the same distance from every vertex 
of a n . When a n is embedded in M. m for m > n, c(a n ) is the unique point that among all 
points equidistant from the vertices of a n minimizes the distance to the vertices of <r™. The 
circumradius of a simplex a, i.e., the distance from c(a) to the vertices of a, is denoted R(a). 

We also use the cone operation of algebraic topology [9], writing u*a n to indicate the simplex 
formed by taking the convex hull of a vertex u together with the n-dimcnsional simplex o~ n to 
form a simplex of dimension n + 1 . This notation may also be used for a set K of simplices; 
u * K is the set of simplices {u* a : a E K}. The affine hull of a set S C R m , which we denote 
aff(5), is the smallest affine space that contains S. For a simplex a n = [v . . . v n ] we can define 
it as 



The affine hull of a simplex a may also be called the plane of a. 

When referring to a simplex er, the boundary of the simplex, denoted Bd(a), is the union of 
the set of proper faces of a, i.e., the set of all faces of a other than a itself. The interior of the 
simplex, denoted by Int(cr), is defined as a — Bd(a). More generally, for a set S, we use Int(S*) 
to refer to the interior of £ taken with respect to the usual topology of aff(S'). For the closure 
of a set S we use the notation Cl(S'). 

For a vertex u of a simplicial complex we define St u, the star of the vertex, to be the union 
of the interiors of all the simplices for which u is a vertex. The closure of the star, or the closed 
star, Cl(St u), is the union of all simplices incident to u. The link of a vertex u is defined by 
Lk u = Cl(St u) — St u. Many of the terms briefly defined here are defined and discussed at 
more length in [9]. 

We wish to avoid any ambiguity about the dimension of the circumsphere or circumball of 
a simplex a. Throughout this paper the objects circumsphere and circumball always live in 
aff(cr). Thus the circumsphere of a triangle is always a copy of S 1 , even when the triangle is 
embedded in R 3 as the facet of a tetrahedron. The equatorial ball of a simplex a, sometimes 
denoted B(a), is a ball of radius R(a) centered at c(cr), but distinguished from the circumball of 
a by the fact that the equatorial ball is considered in a higher-dimensional space. For example, 
the equatorial ball of a triangle r considered in M 3 is the unique 3-dimensional ball that has the 
circumcircle of t as an equator (see Fig. 1). The ambient higher-dimensional space that contains 
the equatorial ball should be made clear wherever the term is used. In this paper, wherever the 
term equatorial ball is used and a higher-dimensional space is not explicitly specified, the term 
appears in the context of a simplex a n = u * r n_1 , and _B(r™ _1 ) is an n-dimcnsional subset of 



We frequently discuss the facets of a simplex a n = [vq . . . v n ] . As a matter of convention, 
the facets usually are denoted t^ -1 , . . . , t" -1 with the understanding that the facet r™ _1 is the 
facet opposite vertex w^. Thus a n — Vi * t™ -1 for each i = 0, 1, . . . , n. 

4 Characterizing the Well- Centered Simplex 

We now investigate geometric properties of an n-wcll-centered n-simplex. The context for this 
discussion is a simplex a n — u * r" -1 with facet r" -1 given. The vertex u is free to move, and 
we wish to determine whether a n is n-well-centered based on the position of u relative to the 
fixed vertices of r n_1 . 




aff(a"). 
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Figure 1: Four views of the same 3-well-centered tetrahedron a 3 in the same orientation. 
From left to right the views show a 3 with the equatorial balls of its bottom, right, left, and rear 
facets. For each facet rf of a 3 , the circumsphere (i.e., circumcircle) of rf , which is an equator 
of the equatorial ball B(rf), is shown. Because the tetrahedron is 3-well-centered, the vertex 
Vi opposite facet t? lies outside of B(rf); an n-simplex is n- well-centered if and only if for each 
vertex vi, Vi lies outside of the equatorial ball of the facet t™ -1 opposite i;,. For the bottom 
facet and rear facet views in the figure, the reader may need to look closely to see that the edges 
incident to Vi do pierce B(rf), and Vi lies outside B(rf). 

We first recall an alternate geometric characterization of the n-well-centered n-simplex and 
state its consequences in this context. The alternate characterization is stated in terms of 
equatorial balls. Using the result, which adopts the notational convention that a n = Vi * r™ _1 
for each i — 0,1,..., n, one can determine whether an n-simplex is n-well-centered without 
explicitly computing c(a n ). 

Theorem 1 (Equatorial Balls Condition). The simplex a n is n-well-centered if and only if 
vertex Vi lies strictly outside ^(r™ -1 ) for each i = 0, 1, . . . , n. 

Proof. See [17]. □ 

Figure 1 illustrates Theorem 1 as it applies to a tetrahedron. For each vertex Vi of the 
tetrahedron, Fig. 1 shows the equatorial ball B(n) of the facet Tj opposite Vi, emphasizing in 
a darker color the circumcircle of r, (which is an equator of B(n)). The figure shows that in 
each case Vi is outside the equatorial ball of Tj, so we can conclude that the tetrahedron is 
3-well-centered. Moreover, this same condition is satisfied by every 3-well-centered tetrahedron. 

In the context of an n-simplex cr™ with a free vertex u opposite a fixed facet r" _1 , Theorem 1 
becomes a necessary condition that vertex u must satisfy if cr™ is to be n- well-centered. 

Corollary 2 (One-Facet Equatorial Ball Condition). Let cr™ = u * t™ -1 . If the simplex a n is 
n-well-centered, then u lies strictly outside o/i?(T n ~ 1 ). 

To introduce the remaining results of this section and get a somewhat different perspective 
on Corollary 2, we consider the sketch in Fig. 2. In the sketch, t is a given triangle in M 3 with 
fixed vertices. Triangle r represents a facet of a tetrahedron a whose fourth vertex u has not yet 
been determined. (Thus tetrahedron a docs not appear in Fig. 2.) We suppose, however, that 
the circumcenter c(a) of the tetrahedron is known, so u is constrained to lie on the circumsphere 
of a. The two sides of Fig. 2 show two different cases; on the left r is not 2- well-centered, and 
on the right r is 2-well-centered. 

Now t lies in the plane aff(r), and the reflections of r and aff(r) through c(a) are t' and 
aff (t') respectively. The plane aff (r) intersects the circumsphere to determine a lower spherical 
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Figure 2: A facet t of a tetrahedron a (not shown in the figure) and the reflection t' of r 
through c(a). In the case on the left r is not 2- well-centered; on the right r is 2-well-centered. 
The tetrahedron a is 3-wcll-centered if and only if its vertex u opposite t lies in the spherical 
triangle determined by the intersection of the circumsphere of a and a cone on t' with apex 
c(a). (The spherical triangle is shown in the case on the right, but not on the left.) 



cup C, and aff (t') determines an upper spherical cup C . The necessary condition of Corollary 2 
says that when a is 3- well-centered u does not lie in the lower spherical cup C. From the geometry 
one can see that, in fact, if a is to be 3- well-centered, then u must lie strictly inside the spherical 
triangle determined by the intersection of C", the upper cup of the circumsphere of a, with the 
geometric cone on t' with apex c(tr). (This spherical triangle is drawn in the case on the right, 
but not in the case on the left.) 

In particular, there is a necessary condition that u must lie strictly in the upper cup C 
of the sphere. Thus, speaking with regard to the orthogonal projection into aff (r), u projects 
(vertically in the figure) to the interior of the circumdisk of r. Moreover, if r is 2-well-centered, as 
on the right, then the projection of the spherical triangle into aff (t) contains the projection of r' 
into aff (r). Thus if the projection of u into aff (t) lies inside the projection of r' into aff(r), this 
is sufficient to establish that a is 3- well-centered. These conditions and their generalizations into 
higher dimensions are the first two conditions discussed in this section. The geometric intuition 
developed here is formalized and proved algebraically in Propositions 3 and 8. 

Finally, we consider varying the position of c(a). Notice that as c(a) moves in Fig. 2 from 
the circumcenter of t upward along a line orthogonal to aff(r), the spherical triangle of it- 
positions that produce a 3-well-centered tetrahedron with circumcenter c(<r) sweeps out a solid 
3-dimensional region. Tetrahedron a will be 3-well-centered if and only if u lies in this region. 
The section closes by describing this region for arbitrary dimensions in terms of polynomial 
inequalities. (See Fig. 9.) 

Now we state a proposition that gives a necessary condition for an n-simplex a n to be 
n- well-centered. See Fig. 3. 



Proposition 3 (Cylinder Condition). Let a n be an n-well- centered n-simplex in M. n with u a 
vertex of a n and r™ _1 the facet of a n opposite u. That is, let a n = u * t™ _1 . Let P be the 
orthogonal projection P : W 1 — > aff(r™ _1 ). Then \\P{u) — c(r™ _1 )|| < i?(r™ _1 ), i.e., u projects 
to the interior of the circumball of r n ~ 1 . 

Proof. Consider the coordinate system on K" such that c{a n ) is the origin, and aff(r" _1 ) = 
{x e R" : x n = k} for some constant k < 0. In this coordinate system, P is the projection 

P : (x\, . . . , x n -i, x n ) i > (x\, . . . , x n _i, k). 
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Let u — (x\, . . . , x n ) in this coordinate system. We have assumed that ct" is n-well-centered, so 
c(a n ) (the origin) is strictly interior to ct". It follows that k < and x n > 0. 

Consider the line segment £ from u to — u. Observe that £ is a diameter of the circumsphere 
of ct". Moreover, t n Int(r™ _1 ) ^ 0. This follows from the fact that ct" is n-well-centered; we 
have ct" = u* r" -1 and c(a n ) G Int(<7™), so there must be some point w € Int(r n_1 ) such that 
c(ct") lies on uio C £ We notice, then, that the point —it lies below aff(Y™ -1 ) and conclude that 
x n ^> k . 

By the Pythagorean theorem, _R(r n_1 ) 2 + k 2 = R(a n ) 2 . We also have 

n 

»=i 

since u lies on the circumsphere of ct". It follows that 

ra-l 

||P(«) - cCr- 1 )!! 2 = J2 x i= R ^f ~ x l< R (° n ? -k 2 = Rir"- 1 ) 2 . 

i=l 

□ 



The statement is not limited to ct" C M", of course. For a simplex ct" C K m with m > n, 
there exists a coordinate system such that 

aff(cr") ={iet m :i i = 0for! = n+l,...,m}, 

and the same proof applies. 

Remark 4. Given a particular simplex t" _1 C R n , Proposition 3 provides a geometric necessary 
condition on the location of vertex u to create an n-well-ccntered simplex ct" = u* T n ~ 1 . Vertex 
u must lie within a solid right spherical cylinder over the circumsphere of r" _1 if ct" is to 
be n-well-centered. Figure 4 illustrates the condition in 2D and 3D, making it clear how this 
condition generalizes from the familiar 2-D case into higher dimensions. In each case the vertices 
of the base simplex r" _1 , as well as the circumcenter c(r" _1 ), are marked by small dark-colored 
balls. If u * T n ~ x is n- well-centered, then the vertex u must lie inside the gray cylinder over 
the circumsphere of r" _1 . In the notation of Fig. 2, where the circumcenter of ct is known, the 
Cylinder Condition says that u must lie either in the upper cup C or the lower cup C . 

Remark 5. The One-Facet Equatorial Ball Condition (Corollary 2) is also a necessary condition 
on the location of vertex u. In M 2 , the combination of Corollary 2 and the Cylinder Condition 
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Figure 4: If a simplex u * r" 1 is n- well-centered, then u is interior to a solid right spherical 
cylinder over the circumsphere of r™" 1 . 

is sufficient to guarantee that a triangle (a 2-simplex) is acute (is 2- well-centered). In M. n for 
n > 3, however, an n-simplex u * r™" 1 for which u satisfies both of these necessary conditions 
might not be n- well-centered. 

Example 6. For example, consider the tetrahedron a = a 3 with vertices (—0.152, 0.864, —0.48), 
(-0.64, -0.6, -0.48), (0.6, -0.64, -0.48), and (-0.192, -0.64, 0.744), whose circumcenter lies at 
the origin. For three of the four facets rf of a 3 , vertex Vi satisfies both necessary conditions with 
respect to jf, and for the fourth facet Vi satisfies the Cylinder Condition, but not the One-Facet 
Equatorial Ball Condition. The tetrahedron a 3 is not 3- well-centered. 

Figure 5 shows several different views of a 3 . The large view at left shows that a 3 is not 3- well- 
centered; the circumcenter of a 3 , marked by a small axes indicator, lies outside the tetrahedron. 
The four small views on the right side of Fig. 5 are views directly down onto the facets rf of 
a 3 . In each case the circumcircle of rf is rendered in a darker color, and one can see that the 
vertex above the facet projects to the interior of the circumball of the facet, i.e., that vertex 
Vi satisfies the Cylinder Condition with respect to t? . In three of the four cases — all except 
the case at lower left — the vertex also satisfies the One-Facet Equatorial Ball Condition. The 
particular example in Fig. 5 is also mentioned in [15], which gives some additional statistics on 
the tetrahedron. 

Example 7. The tetrahedron with vertices at (-0.01,-0.01,-0.01), (1,0,0), (0,1,0), and 
(0,0, 1) is another tetrahedron that is not 3- well-centered. It also has three vertices that satisfy 
the One-Facet Equatorial Ball Condition and four vertices that satisfy the Cylinder Condition. 
This example is dihedral acute, in contrast to the previous example. 

The above examples illustrate that the One-Facet Equatorial Ball Condition and the Cylinder 
Condition are not enough to establish that the n-simplex u * r' i_1 is n- well-centered. However, 
the following proposition does provide sufficient conditions that u* r™" 1 is n- well-centered. See 
also Fig. 6. 

Proposition 8 (Prism Condition). Let T n ~ x be an (n — 1) -well-centered simplex in R n and 
<r n = u * T n ~ . If u lies outside the equatorial ball B(r" _1 ) and the reflection of P(u) through 
c(t™ -1 ) j s interior to T n ~ x , then a n is n-well- centered. 

Proof. We assume the stated hypothesis and take the same coordinate system that was used in 
the proof of Proposition 3. Observe that if u were on the equatorial ball of r™" 1 , then c(a n ) 
would lie in r" _1 , coinciding with c(r"~ 1 ). Because u lies outside the equatorial ball of r™ _1 , 
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Figure 5: A tetrahedron that is not 3- well-centered, even though every vertex satisfies the necessary 
condition of Proposition 3. Three of the vertices also lie outside the equatorial balls of their respective 
opposite facets. 



Figure 6: Because P(—u) and c(r" 1 ) = P(c(a n )) are both interior to r™ 1 and c(a n ) is above 
aff (r n_1 ), we know that the tetrahedron a n = u * t 11-1 is 3- well-centered. 

c(<7 n ) lies interior to the same halfspace as u with respect to aff(r™ _1 ). It follows that k < 
and x n > k. 

Observe that, as shown in Fig. 6, the reflection of P(u) through c(r n_1 ) is P(—u). By the 
hypothesis, P(—u) is interior to T n . Thus P{—u) is interior to the circumball of a n and 



It follows that \x n \ > \k\ = —k. Since we know that x n > k, we conclude that x n > — k > 0. 

Let I be the line segment from u to —u. We will show that I intersects the interior of T n_1 . 
Then, because a n = u * r n ~ 1 and k < < x n (so that c(a n ) £ I is above r n_1 and below u), 
we will be able to conclude that c(a n ) is interior to er™. We know that P(c(a n )) = c(r n_1 ) 
is interior to t™ _1 because r n_1 is (n — l)-well-centered. Since P(—u) and P(c(a n )) are both 
interior to r™ _1 , the line segment from c(er n ) to —u, which is contained in t, is interior to the 
(convex) infinite prism r n_1 x R. Moreover, > k > —x n (i.e., c(<7™) is above r™ _1 and — u is 
below r n_1 ), so this part of segment £ intersects the interior of r™ _1 . □ 




n-l 



n 



\\p(u)\\ 2 = \\p(-u)\\ 2 = e + Y, x *< i? ( fT ") 2 = E x * 



i=i i=i 



8 



Figure 7: A tetrahedron for which the top vertex and bottom facet satisfy all of the sufficient 
conditions for being 3-well-centered except that the bottom facet is not 2- well-centered. 

As was the case for Proposition 3, Proposition 8 is not limited to a n C R"; in higher-dimen- 
sional spaces K m there is a coordinate system such that 

aff (cr n ) = {x e W n : Xi = for i = n + 1, . . . , m}, 

and the same proof applies. 

After reading Proposition 8, one might ask whether the requirement that the facet r n_1 be 
(n — l)-well-centered can be removed from the proposition. It may already be clear from the 
discussion of Fig. 2 that the answer to this question is no. The tetrahedron in Fig. 7 is an 
explicit example that confirms the requirement cannot be removed. 

Example 9. The tetrahedron in Fig. 7 is the convex hull of vertices Vq — (0.224, —0.768, —0.6), 
Vi = (0.8,0,-0.6), v 2 = (0.224,0.768,-0.6), and v 3 = (-0.28,0,0.96). The bottom facet in 
Fig. 7, which is the triangle r| — [VQV1V2}, lies in the plane £3 = —0.6 and is an obtuse triangle. 
The obtuse angle is at vertex V2, the rightmost vertex in Fig. 7. Taking this bottom facet to 
be t 2 as in Proposition 8, and the top vertex to be u — 1*3, we satisfy the conditions that u 
lie outside the equatorial ball of r 2 and that the reflection of P(u) through c(t 2 ) be interior 
to r 2 . Indeed, c(r 2 ) = (0,0,-0.6) and i?(r 2 ) = 0.8 with \\u - c(r 2 ) || = V2ll2 w 1.58, so u 
is outside B(t 2 ), and P(—u) = (0.28,0,-0.6) lies inside r 2 . Thus we satisfy all of the Prism 
Condition except the requirement that r 2 be 2- well-centered. It is clear from Fig. 7 that this 
is not sufficient; the circumcenter of tetrahedron u * r 2 , marked by a small axes indicator, lies 
outside the tetrahedron, so the tetrahedron is not 3-well-centered. 

Remark 10. Like the condition of Proposition 3, the condition of Proposition 8 has a nice 
geometric interpretation. Given an (n — l)-well-centered facet r n_1 , if the vertex u opposite 
t™ -1 lies outside i?(r"~ 1 ) and within an infinite prism (a right cylinder) over the reflection of 
r™ -1 through its circumcenter, then <r n — u * t™ _1 is n- well-centered. Figure 8 portrays the 
region defined by the Prism Condition for specific examples in 2 and 3 dimensions. In each case 
the base simplex r™ -1 is shown in dark colors and solid lines, and its reflection is outlined with 
lighter colors and dashed lines. In the figure, each r™ _1 is (n — l)-well-centered, so for a vertex 
u lying inside the prism over the reflection of t" _1 through c(t™ -1 ) and outside the equatorial 
ball of r" _1 , i.e., for a vertex u lying in the gray region shown in Fig. 8, the simplex u * r™ _1 
will be n-well-centered. Note that on the left in Fig. 8 the base simplex and its reflection should 
actually lie on top of each other, but are set slightly apart in the drawing so the reader can 
distinguish them from each other. 

We have now established two different conditions for an n-simplex to be n-well-centered. 
One condition is a necessary condition, and the other condition is a sufficient condition. Both 
conditions are stated in terms of the location of a vertex u relative to the facet t"" 1 opposite u. 
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Figure 8: If the base simplex r" -1 is (n — l)-well-centered and vertex u is both outside i?(r n_1 ) 
and inside the infinite prism over the reflection of r™ _1 through c(r™~ 1 ), then the simplex u * r n_1 is 
n- well-centered. 

The regions defined by the necessary condition and the sufficient condition may be quite different 
from each other. For example, in the 3-D portions of Figs. 4 and 8 the same base simplex r™ _1 
yields rather different regions for the two conditions. It is natural to seek a precise description of 
the region where the vertex u will produce an n- well-centered n-simplex u*r" _1 . The following 
discussion develops just such a set of conditions on the location of u. The conditions take the 
form of a system of cubic polynomial inequalities in the coordinates of u. The simplex u * r™ _1 
will be n- well-centered if and only if the coordinates of u satisfy the polynomial inequalities. 

The inequalities are derived from a linear system of equations discussed in [1]. This linear 
system, which provides one way to compute the circumcenter of a simplex o~ n embedded in W n for 
m > n, is briefly reviewed here. We may write the circumcenter c of a simplex a n = [vqv-l . . . v n ] 
as a linear combination of the vertices Vi £ M. m , 



c = a a v + a 1 v 1 H h a n v n , 

with the coefficients at satisfying X)™=o ai = ^ ne coefficients on are known as the bary 'cen- 
tric coordinates of the circumcenter. The condition that a n be n-well-centered is the same as 
the condition that < aj for every a,, i.e., the condition that the circumcenter be a convex 
combination of the vertices of a n with strictly positive coefficients. 

Suppose we are given the coordinates of the vertices Vi of a n . We know that 

(c — Vi,c - Vi) = ||c - wj 2 = R 2 

for each vertex Wj. Introducing the variable A = R 2 — ||c|| 2 , we obtain the n + 1 equations 
2(c, Vf) + A = ||y;|| 2 . Since the vertices Vi are known, each equation is a linear equation in 
the n + 2 unknowns oto, «i, . . . , a n , A. The final equation of the system is X)T=o ai ~ ■"•> wmc ^ 1 
forces the ai to be barycentric coordinates. As long as this linear system of n + 2 equations in 
n + 2 unknowns is nonsingular, we can solve for the barycentric coordinates. If the simplex is 
nondegenerate, i.e., if the n + 1 vertices are affinely independent, then the simplex has a unique, 
finite circumcenter, which has unique barycentric coordinates. It follows that the linear system 
has a unique solution; hence the matrix is nonsingular. 
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Let A be the matrix of this linear system and b the right-hand side, 



A = 



(2{v ,v ) 2{v ,vij 

2(v 1 ,v ) 2{vi,vi) 

2(v n ,v ) 2(v n ,vi) 

V i i 



2(v ,v n ) 1\ 

2(ui,U n ) 1 

2(v n ,v n ) 1 

1 0/ 



b = 



( (vo,v )\ 
{v n ,v n ) 

V 1 / 



For i = 0, 1, . . . ,n we let be the matrix A with column i + 1 replaced by 6. Cramer's rule 
tells us that on = det(Aj)/det(A). If we consider vertices vo, ■ ■ ■ ,v n -\ to be the vertices of 
some given r™ _1 and v n to be a free vertex it, then the barycentric coordinates on are rational 
functions of the coordinates of u. Thus the conditions on > become algebraic inequalities in 
the coordinates of u. 

To simplify matrix A a little, we translate each vertex of the simplex by —vq- The translation 
may change the value of A in the solution vector in fact, A = always holds for the translated 
system — but the barycentric coordinates of the circumcenter are not changed by translating 
the vertices of the simplex. If m > n we make one further simplification. In the translated 
coordinate system, we rotate the simplex about the origin v to obtain a simplex for which 
vector Vi — vq G {x : x n +\ = ■ ■ ■ = x m = 0} for each i = 1, . . . , n. Rotation about the origin is 
an orthogonal transformation, so it does not change any of the entries of the linear system and 
does not affect the barycentric coordinates. 

If we restrict our attention to one of the open halfspaces bounded by aff(r™ _1 ), we have 
either det(A) > or det(A) < throughout the halfspace, because det(A) is a continuous 
function of the entries in A and A is singular only when u G aff(T n_1 ). We will see that, in fact, 
det(A) < holds everywhere, so det(A) < throughout the halfspace. 

The first row and the first column of A in the simplified linear system are all zeroes except 
for the last entry, which is 1 in both cases. Computing the determinant of A by first expanding 
across the first row and then expanding down the first column (one with an odd number of 
entries and the other with an even number of entries) we find that det(A) = — det(_B) where B 
is the submatrix of A spanning rows 2 to n + 1 and columns 2 to n + 1. The n x n submatrix 
B has the form 2V T V, where V is the m x n matrix 



V = (yi - v v 2 - v 



vo) ■ 



Because of the earlier rotation of the simplex, the last m~ n coordinates of each vector Vi — vo 
are zeroes, and if we take V to be the first -n rows of V, then V is an n x n matrix that satisfies 
V T V^= V T V. It follows that B = 2V T V. Thus dct(B) = 2™det(y) 2 > 0. Observing that 
det(y) is the signed volume of the parallelepiped spanned by the vectors that form the columns 
of V, we note that det(-B) > holds when the columns of V are linearly independent, i.e. when 
the vertices of the original simplex arc affincly independent. 

Thus with the assumption that r" _1 is a fully (n — l)-dimensional simplex, we know that 
det(A) < when the vertex u lies in either of the open halfspaces bounded by aff (r™ _1 ). For u 
outside aff(T n_1 ), then, we conclude that oti = det(Aj)/ det(A) > if and only if det(Aj) < 0. 
Hence the simplex u * r™ _1 will be n-well-centered if and only if the coordinates of u satisfy the 
polynomial inequality det(-Aj) < 0. 

It remains to show that the equation det(Aj) = is a polynomial in the coordinates of u 
of degree at most 3. To do this we examine the entries of Ai that depend on u. All of these 
entries appear in row n + 1 or in column n + 1. At most two of these entries are quadratic in 
the coordinates of u — the entry at position (n + 1, n + 1) and the entry at (n + 1, i + 1). (Only 
one entry is quadratic in the coordinates of u when i = n.) Every other entry that depends on 
u is linear in the coordinates of u. Using S n to denote the group of permutations on n letters, 
the determinant of an n x n matrix M can be written as 



det(M) = J2 s S n ( 7r ) II 



(:))• 



7TSS„ 
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Figure 9: Given a facet r n_1 , the region where the vertex u may lie to produce an n-well- 
centered simplex u * r™ _1 is defined by a system of polynomial inequalities. When r" _1 is 
(n — l)-well-centered, so that there are regions related to both the necessary and sufficient 
conditions, the actual region where u may lie is somewhere in between the regions defined by 
the necessary Cylinder Condition (Fig. 4) and the sufficient Prism Condition (Fig. 8). 



where Mjk stands for the entry in row j, column k of matrix M, and sgn(7r) is the signum 
function applied to the permutation. Considering the structure of matrix Ai, we observe that 
each product in this definition of det(^li) involves at most two terms that depend on u, and 
at most one of these — the entry selected from row n + 1 — is quadratic in the coordinates of u. 
Thus the determinant is a summation of terms that are polynomial in the coordinates of u and 
have degree at most 3. 

We can also explain this from the perspective of computing the determinant by expanding it 
along a row or column. We will consider a specific example with i = 2 arising from a tetrahedron 
(dimension n — 3), but the discussion applies to the general case. We have 



(2{v ,Vq) 

2(wi,uo) 

2(V2,K ) 

2(u, v ) 

V i 



2(«o,fi) 

2(»i,T>i) 

2<«,«i) 
1 



(vo,vq) 

(Vl,Vi) 
(V2,V 2 ) 

(u,u) 

1 



2(v ,u) 
2{ Vl ,u) 
2{v 2 ,u) 
2{u,u) 
1 



1\ 
1 
1 
1 

V 



for this particular example. If we compute det(Ai) by expanding down column n+l (column 4, 
in this case), we find that term n + l of the summation is a quadratic function of the coordinates 
of u multiplied by the determinant of a submatrix that is constant with respect to u. In our 
example, this is the fourth term in the summation, 



2(u, u) ■ det 



/2{v ,v ) 
2(«i,fo) 

\ 1 



2(«o,«i) 
2{v 1 ,v 1 ) 
2{v 2 ,v 1 ) 
1 



(v ,v ) 1^ 

(vi,vi) 1 

(V2,V 2 ) 1 

1 0/ 



The remaining n + 1 of the terms in the summation are linear (or constant) functions of u 
multiplied by a determinant of some other submatrix of Ai that is not constant with respect to 
u. For our example, the first term of the summation is 



-2(v 0l u) ■ det 



(2{vx,v Q ) 2{vi,v x ) 

2(v 2 ,v ) 2(v 2 ,vi) 

2(u,v ) 2(u,vi) 

V i i 



(Vi,Vi) l\ 

(v 2 ,v 2 ) 1 

(u,u) 1 

i oy 
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Expanding the appropriate row (usually row n) of each of these submatrices in similar fashion, 
we obtain a summation of terms that are either linear or quadratic in the coordinates of u (at 
most one term is quadratic) , each multiplied by the determinant of a smaller submatrix that is 
constant with respect to u. 

We state the conclusions of the foregoing discussion as a formal proposition. 

Proposition 11. Let <r n = m*t" _1 for a fixed facet t" -1 . The n-simplex a n is n-well-centered 
if and only if the coordinates of vertex u satisfy the inequalities det(-Aj) < 0, which are cubic 
polynomial inequalities in the coordinates of u. 

Figure 9 gives a graphical representation of the precise region where the vertex u may be 
placed to produce a 3- well-centered tetrahedron u * t" _1 . The facet t"^ 1 used in Fig. 9 is the 
same facet used to illustrate the necessary condition for a tetrahedron in Fig. 4 and the sufficient 
condition for a tetrahedron in Fig. 8, so readers can see for this specific case how the full region 
compares to the regions defined by the necessary condition and the sufficient condition. The 
facet t™ -1 along with its circumcircle and the reflection of r™ _1 through c(t™ _1 ) are shown in 
Fig. 9 to aid this comparison. It should also be noted that Fig. 9 was generated using MATLAB's 
isosurfacc function and evaluations of the polynomial inequalities on a finite grid, so the graphical 
representation has some slight imperfections. For instance, the entire circumcircle of r™ _1 lies 
in the boundary of the region even though in Fig. 9 it appears that there is a small gap above 
and below aff^™" 1 ). 

5 Local Combinatorial Properties of 3- Well-Centered Tetra- 
hedral Meshes 

The geometric properties of the n-well-centered n-simplex discussed in Sec. 4 have implications 
for the combinatorial properties of well-centered meshes. As a simple motivating example we 
consider the 2-dimensional case of a triangle mesh in the plane. If v is a vertex interior to this 
mesh and there are fewer than five edges incident to v, then some angle incident to v has measure 
7r/2 radians or larger. Thus the mesh has a nonacute triangle. This geometric observation can 
be restated as a combinatorial property of 2- well-centered (i.e., acute) triangle meshes. Namely, 
there are at least five edges incident to every interior vertex of an acute triangle mesh in R 2 . 
This well-known fact is a key ingredient in the generation of 2-wcll-ccntered triangle meshes 
through optimization of vertex coordinates; the mesh must satisfy this combinatorial condition 
at every interior vertex if optimizing the vertex coordinates is to have any hope of finding an 
acute mesh. 

Similarly, tetrahedral meshes in R 3 that are 2-well-centered or 3-well-centered must satisfy 
certain local mesh connectivity conditions. These combinatorial conditions, which are key to 
creating well-centered tetrahedral meshes, are analyzed in the next two sections. This section 
develops some of the combinatorial properties of 3-well-centered tetrahedral meshes, and the 
next section examines combinatorial properties of 2-well-centered tetrahedral meshes. 

The combinatorial properties of tetrahedral meshes in R 3 are more complex than the anal- 
ogous properties for triangle meshes in R 2 . In a triangle mesh in R 2 , the link of an interior 
vertex is a set of edges that form a cycle around the vertex, i.e., a triangulation of a topolog- 
ical circle (5* 1 ). The number of edges incident to the interior vertex, which is the number of 
vertices on the cycle, completely characterizes the neighborhood of the vertex. In tetrahedral 
meshes in R 3 , on the other hand, the link of an interior vertex is a triangulation of a topological 
sphere S 2 . Thus the number of edges incident to the vertex does not completely characterize 
the neighborhood of the vertex. We do, however, prove necessary conditions on the number 
of edges that must be incident to an interior vertex in a tetrahedral mesh in R 3 in order for 
the mesh to be 3-well-centered, 2-well-centered, or completely well-centered. We also show that 
there is no sufficient condition in terms of the number of edges incident to an interior vertex. 
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Figure 10: If Vi and Vj are both contained in the interior of the solid right spherical cylinder over 
the circumcircle of t u , then u lies outside the simplex formed from m, Vj, and the other vertices of t u . 



Much of the discussion in Sees. 5 and 6, then, is phrased in terms of the link of an interior 
vertex. For a tetrahedral mesh in M 3 , this is a triangulation of S 2 , which corresponds to a 
planar triangulation in a graph theoretic sense. We try to avoid the term planar triangulation 
to prevent possible confusion with triangle meshes in R 2 . 

We begin with two results that apply in arbitrary dimension. The first lemma generalizes the 
following statement about planar triangle meshes, using the Cylinder Condition (Proposition 3) 
to relate geometry to combinatorics. If a planar triangle a 2 = [abc] is subdivided into three 
triangles by adding a vertex u interior to a 2 and adding edges [ua], [ub], and [uc] to obtain 
u * Bd(cr 2 ), then at most one of the three triangles [abu], [bcu], and [cau] in u * Bd(cr 2 ) is an 
acute triangle. The main idea of the proof of Lemma 12 is illustrated by Fig. 10. 

Lemma 12. For n > 2, let a n = [vqVi . . . v n ] have facets Tq~ , r™ _ , . . . , t™ _1 . If u is a point 
lying in a n , then at most one of the n-simplices of u*Bd(a n ), i.e., at most one of the simplices 
u * Tq~ X ,u * t" -1 , . . . ,u* t™ -1 , is an n-well- centered n-simplex. 

Proof. It suffices to prove the statement when u is in the interior of a n . Indeed, if u is on the 
boundary of o~ n and two or more of the simplices u * r™" 1 are n- well-centered, then we can 
slightly perturb u into the interior and obtain a point u e Int(<7™) with at least two n-well- 
centered n-simplices. Thus we assume that u € Int(cr n ). 

Let t™ -1 and t™ -1 be two distinct facets of u n . Then u * t™ _1 and u * t"™" 1 are n-simplices, 
and r™ -1 n t™ -1 is an (n - 2)-dimcnsional face of a n . The face t" -1 n r"" 1 is incident to all 
but two of the vertices of cr", the two vertices Vi and Vj. (Recall that Vi is opposite r™ _1 and Vj 
is opposite t™ -1 .) Notice that u * r™ _1 and u * t™ -1 have a common facet, the (n — l)-simplex 
r™" 1 := u * (r/ l_1 n Wc let T C aff(er™) be the solid right spherical cylinder over the 

circumball of r™ _1 . 

Assume towards contradiction that u * r™ _1 and u * are both n-well-centered. By the 
Cylinder Condition (Proposition 3), both and Vj lie in Int(T). Now T is a convex set, and all 
the vertices of a n lie in T, so a n C T. On the other hand, u lies on the circumsphcrc of r™ _1 , 
so u e Bd(T). Thus u Int(cr") C Int(T), contradicting the assumption we made in the first 
paragraph of the proof. We conclude that at most one of u*r™ _1 , u*r" _1 is n-well-centered. □ 



The next theorem shows that Lemma 12 has implications for the local combinatorial proper- 
ties of n-wcll-ccntcrcd meshes. The theorem is stated using the language of simplicial complexes. 
We say that a vertex u is an interior vertex in an n-dimcnsional simplicial complex embedded 
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Abstract Complex K, Embedding of Induced Geometric 

Bd(K) = L M = u*L Realization of K 



Figure 11: In R 2 the simplicial complex K consisting of two triangles and [i>;lI>31>2] and their 

faces (left) satisfies the hypothesis of Theorem 13 for the mesh M — u* L embedded in R 2 (center), 
so M is not 2- well-centered. The embedding of M in R 2 induces a geometric realization of K into 
R 2 (right). The geometric realization of K is not an embedding in this case, since [K1U3U2] is inverted 
here. The particular embedding of M does not affect the existence of the abstract complex K in 
Theorem 13; there is no embedding of M that is 2-well-centered. 



in R n if (the underlying space of) Lk u is homeomorphic to 5 n_1 , the sphere of dimension n—1. 
Thus the closed star of u is homeomorphic to an n-dimensional ball in R n , and the point v lies 
in the interior of the ball in the standard topology on R™. 

When we speak of an abstract simplicial complex K we make an important distinction 
between an embedding of K and a geometric realization of K. An embedding of K is an 
assignment of coordinates in R™ to the vertices of K such that K is a simplicial complex in R™ 
with vertices at the specified locations. By a geometric realization of K we mean merely some 
assignment of coordinates in R™ to the vertices of K. Thus in a geometric realization of K in 
R™, it is possible for K to have self-intersections. Figure 11, which is related to the proof of 
Theorem 13, illustrates the distinction between these two terms. 

Theorem 13 (One-Ring Necessary Condition). Let u be an interior vertex of an n-dimensional 
simplicial complex M (e.g., a mesh) embedded in R™, and set L = Lk u. If there exists an 
abstract finite n-dimensional simplicial complex K such that 

(i) K is an n-manifold complex (with boundary) 

(ii) Bd(K) is isomorphic to L, and 

(Hi) for every n-simplex a n G K, there are at least two (n — l)-simplices in Bd(a n ) n L, 
then u* L is not n-well-centered. 

Proof. We first observe that every vertex of K must also be a vertex of L. By assumption (i), 
every simplex of if is a face of some n-simplex of K, so if K had a vertex v not in L, then 
there would be some n-dimensional simplex a n € K incident to v, and a n would have only 
one (n — l)-dimcnsional face not incident to v. Since v ^ L, it follows that Bd(er n ) n L would 
contain at most one (n — l)-simplex, and (iii) would not be satisfied. 

The embedding of M in M. n includes an embedding of u * L in R™ . Since every vertex of K 
is a vertex of L, this embedding of u * L in R™ induces a geometric realization of K into R™. 
(As shown in Fig. 11, the geometric realization might not be an embedding.) 

We have an embedding of the simplicial complex u * L in R™. Since it is an embedding, 
each n-dimensional simplex is a fully n-dimensional geometric object, and we have consistent 
orientation. Moreover, L is star-shaped with respect to u. We claim that by (i) and (ii) this 
implies that there is some simplex in the induced geometric realization of K that contains the 
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point u (possibly on its boundary). We return to this claim in a moment, but first we show how 
this completes the proof. 

Fix a simplex cr" e K that contains u. Now consider the n-simplices of u * Bd(cr"). By 
assumption (iii) of the hypothesis, at least two of these simplices have a facet in L. Each 
simplex of u * Bd(cr") with a facet in L is a member of u * L, and by Lemma 12 at most one of 
these simplices is n- well-centered. We conclude that at least one of the simplices of u * L C M 
is not n-well-centered. 

Now we prove the claim that there is a simplex of the geometric realization of K that contains 
the point u. Choose a line I through u in general position. General position here means that i 
does not intersect any face of K of dimension less than n — 1. Such an £ can be chosen unless 
u itself lies on a simplex p k of K of dimension k < n — 1, and in that case we are done, since 
there is some cr" >- p k that contains u. 

Since u * L is a simplicial complex and L is star-shaped with respect to u, £ intersects 
exactly two simplices of L, each of dimension n — 1, and the intersection points are in opposite 
directions from u along £. For reference, we designate a + and a — direction and name facet 
t™ -1 (resp. t" _1 ) as the facet of L intersected by I in the + (— ) direction from u. Starting 
from t" -1 we describe a walk along £ through n-simplices and (n — l)-simplices of the geometric 
realization of K that ends at r" _ . By continuity of this walk and r" _ , r™ _1 in opposite 
directions from u, there must be some n-simplex in the geometric realization of K that contains 
u. 

The walk is as follows. Since K is a manifold with a boundary and t™ -1 is on the boundary, 
there is a unique cr™ incident to Tq 1-1 := t™ -1 . Then for a given cr™ the walk is on £ at r™^ 1 , 
and £ intersects some unique second facet of cr™, which we name t™~ . As long as t™ -1 ^ r" _ , 
we are not on the boundary of if, so (since if is manifold) there are exactly two n-dimensional 
simplices incident to r" _1 . One of these is cr™, and the other we name cr™ +1 . Since if is a 
manifold complex, the sequence t™ -1 has no repetitions and must eventually end at r™" 1 . (The 
cr" in the sequence may flip back and forth in orientation, which corresponds to the walk going 
back and forth along £.) □ 



It is worth noting that the existence of the abstract simplicial complex K has no dependence 
on the particular embedding of M in R n . Theorem 13 is really a combinatorial statement, and 
we can use it to show that a particular abstract simplicial complex L — Bd(if ) cannot appear 
as the link of an interior vertex in an n-well-centered mesh embedded in W 1 . 

The case n = 3 is of particular interest. Using the One-Ring Necessary Condition of The- 
orem 13 it is fairly easy to establish a tight lower bound on the number of edges incident to a 
vertex in a 3- well-centered tetrahcdral mesh embedded in R 3 . 

Corollary 14. Let M be a 3 -well- centered tetrahedral mesh embedded in M 3 . For every vertex 
u interior to M, at least 7 edges of M are incident to u. 

Proof. Britton and Dunitz have assembled a catalog of all polyhedra with at most 8 vertices, 
which includes all the triangulations of S 2 with at most 8 vertices [2]. By Theorem 13 it suffices 
to show that each such triangulation L of S* 2 with at most 6 vertices has a corresponding 
tetrahedral complex K such that each tetrahedron of K has at least two facets in common with 
L. 

There is only one triangulation of S 2 with 4 vertices — the boundary of a tetrahedron. The 
corresponding tetrahedral complex is that single tetrahedron. 

There is also only one triangulation of S 2 with 5 vertices. This triangulation is shown in 
Fig. 12 along with two corresponding tetrahedral complexes. Either complex certifies that the 
triangulation cannot be the link of any vertex in a 3-well-centered mesh. 

For six vertices there are two nonisomorphic triangulations of S 2 . The first is shown in 
Fig. 13 along with its corresponding tetrahedral complex. The second is drawn in Fig. 14 along 
with its corresponding tetrahedral complex. □ 
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Figure 12: There is only one triangulation of S 2 with 5 vertices, and it has two corresponding tetra- 
hcdral complexes such that each tetrahedron has at least two facets in common with the triangulation. 





Figure 14: In the other triangulation of S 2 with 6 vertices, the degree list is (5,5,4,4,3,3). This 
triangulation of S 2 also has a corresponding tetrahedral complex such that each tetrahedron has at 
least two facets in common with the triangulation. 
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Figure 15: A 3-well-centered mesh with an interior vertex u such that Lk u has seven vertices 
and degree list (5,5,5,4,4,4,3). The vertex coordinates are listed in the table at right; vertex 
u is at the origin. 
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Figure 16: A 3-well-centered mesh with an interior vertex u such that Lk u has seven vertices 
and degree list (6,5,5,5,3,3,3). The vertex coordinates are listed in the table at right; vertex 
u is at the origin. 



When there are m > 7 vertices, there exist triangulations L of S 2 with m vertices such 
that there is no tetrahedral complex K satisfying both Bd(K) = L and the condition that 
every tetrahedron of K have at least two facets in L. In particular, the triangulations of S 2 
with 7 vertices and degree lists (5,5,5,4,4,4,3) and (6,5,5,5,3,3,3), i.e., polyhedra 7-1 and 

7- 4 in the catalog of Britton and Dunitz, both can appear as the link of a vertex in a 3-well- 
centered mesh. Figure 15 shows an example of a 3-well-centered mesh in M 3 consisting of a 
single vertex u and its neighborhood Cl(St u) such that Lk u is a triangulation with degree 
list (5,5,5,4,4,4,3). Figure 16 shows a similar example for the degree list (6,5,5,5,3,3,3). 
There are three other triangulations of S 2 with 7 vertices. Each has a corresponding tetrahedral 
complex K satisfying the requirements of the One-Ring Necessary Condition (Theorem 13), so 
none of these triangulations can appear as the link of a vertex in a 3-well-centered mesh. 

There are 14 nonisomorphic triangulations of S 2 with 8 vertices. Of these, 5 have tetrahedral 
complexes K that certify they cannot be the link of a vertex in a 3-well-centered tetrahedral mesh 
in K 3 . Each of the other 9 triangulations can appear as the link of a vertex in a 3-well-centered 
tetrahedral mesh in M 3 . (We mention these results without proof here.) For m < 8 vertices, 
then, the necessary condition of Theorem 13 completely characterizes which triangulations can 
and cannot be made 3-well-centered. We leave open the question of whether the One-Ring 
Necessary Condition stated in Theorem 13 is a complete characterization for m > 8 vertices in 
3 dimensions or for n-well-centeredness in dimensions n > 4. 

The triangulations on 8 vertices that cannot be made 3-well-centered are polyhedra 8-4, 8-5, 

8- 6, 8-7, and 8-13 in the catalog [2] of Britton and Dunitz. It is interesting to note that the 
degree list of 8-7, which cannot be made 3-well-centered, is the same as the degree list of 8-8, 
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which can be made 3- well-centered [14]. Thus the degree list of a triangulation does not provide 
enough information to determine whether the triangulation can be the link of a vertex in a 3-well- 
centered tetrahedral mesh in M 3 . There are 50 nonisomorphic triangulations with 9 vertices and 
an exponentially growing number of triangulations with more vertices [8], so although making 
a catalog for 9 or 10 vertices might be somewhat interesting, something more abstract will be 
necessary to definitively characterize which triangulations can be made 3- well-centered. 

In the rest of this section we discuss some more general results in the direction of char- 
acterizing which triangulations of S* 2 can appear as the link of a vertex in a 3-well-centered 
mesh in R 3 . The results fall short of a complete characterization, but do show that the set of 
triangulations of S 2 that cannot appear as the link of a vertex in a 3-well-centered mesh and 
the set of triangulations that can appear as the link of a vertex are both infinite. 

Corollary 15. For any integer to > 4 there is a triangulation of S 2 with to vertices that cannot 
appear as the link of a vertex in a 3-well-centered mesh. 

Proof. We have already proved that this holds for 4 < to < 6. 

For m > 7 we note that the tetrahedral complexes shown on the right hand sides of Figs. 12 
and 13 can be generalized. Consider a tetrahedral complex K consisting of a set of to — 
2 tetrahedra that close around a common edge. The complex K satisfies the conditions of 
Theorem 13, so Bd(if) cannot appear as the link of a vertex in a 3-well-centered mesh. Bd(_ft') 
is a triangulation of S 2 on m vertices with degree list (m — 2, to — 2, 4, . . . , 4). □ 

We note that by removing a single tetrahedron from the example complex K of Corollary 15, 
we obtain another infinite family of triangulations of S 2 that cannot appear as the link of a vertex 
in a 3-well-centered mesh. Each member of this family is a triangulation on to vertices with 
degree list (to — 1, to — 1, 4, 4, 3, 3). This family generalizes the tetrahedral complexes shown 
on the left hand side of Fig. 12 and the right hand side of Fig. 14. 

Corollary 15 is one instance that shows how substantial the difference is between tetrahedral 
and triangle meshes. In the case of triangle meshes in R 2 , where we consider triangulations 
of S 1 as the link of a vertex, the only two triangulations that cannot appear as the link of a 
vertex are the 3-cycle and the 4-cyclc. In contrast, there arc infinitely many triangulations of S 2 
that cannot be the link of a vertex in a 3-well-centered mesh in M 3 . One may wonder whether 
there are still infinitely many triangulations of S 2 that can appear as the link of a vertex in a 
3-well-centered mesh in M 3 . The answer is yes. One way to prove this is to explicitly construct 
an infinite family of 3-well-centered meshes with different vertex links. We will do exactly that 
in a moment, with the help of the following lemma, which we prove using the Prism Condition 
(Proposition 8). 

Lemma 16. Let Sq^ 1 be a unit (n — l)-sphere centered at a point u. If r n ~ l is an (n — 1)- 
well-centered (n — l)-simplex whose vertices lie on iS 1 ™ -1 , and the distance from u to aff(T n_1 ) 
is greater than 1/^/2, then a n := u * r™ _1 is an n-well- centered n-simplex. 

Proof. Suppose that r™ _1 is an (n — l)-well-centered simplex meeting the conditions specified 
in the hypothesis. Let S£ ~ 2 be the circumsphere of t" _1 . Sq ~ 2 is the intersection of aff (t" _1 ) 
with 5q _1 , i.e., an (n — 2)-sphcrc lying in S'q -1 . The orthogonal projection of u into aff(r™ _1 ), 
which we denote by P(u), is the center of Sq~ 2 , i.e., the circumcenter c(r"" 1 ) of t" _1 . 

Since r" _1 is (n — l)-well-centered, it contains the point c(t" _1 ). Thus r™ _1 contains the 
reflection of P(u) through c(r n_1 ). The circumradius of r™ _1 satisfies i?(r™ _1 ) 2 + z 2 = 1, where 
z is the distance from u to aff(T™ _1 ), so because z 2 > 1/2, we have it^r™" 1 ) < 1/V2, and u lies 
outside the equatorial ball of r"" 1 . By the Prism Condition, a n is n- well-centered. □ 
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Figure 17: For fc > 4 we can create an acute triangulation of the unit sphere from a set of 
vertices consisting of the north and south poles and two out-of-phase regular fc-gons. Coning 
such a triangulation to the origin produces a completely well-centered tetrahedral mesh. The 
figure shows the tetrahedral mesh obtained for k = 7. 



It is relatively straightforward to prove the converse as well. For a n — u * r' i_1 with the 
vertices of r n_1 lying on a sphere S'q -1 centered at u, if t 11-1 is not (n — l)-well-centered or the 
distance z from u to aff(r™ _1 ) satisfies z < l/y/2, then cr™ is not n-well-centered. This proof is 
left to the reader; the result is not needed in this paper. 

The simplex a 11 = u * r" _1 in Lemma 16 is an isosceles simplex with all vertices of r ,l_1 
equidistant from the apex vertex u. When n = 2, Lemma 16 reduces to the statement that an 
isosceles triangle is acute if the apex angle is acute. In higher dimensions Lemma 16 tells us 
when an isosceles simplex is Ti-well-centered. Note that in an isosceles simplex all of the faces 
incident to the apex vertex u are isosceles; the plane of each such face intersects the sphere S*™ -1 
in some lower-dimensional sphere centered at u, and Lemma 16 can be applied to these isosceles 
faces. It follows that a n will be completely well-centered if r™ -1 is completely well-centered and 
z > l/v2- lu particular, for the case n = 3, an isosceles tetrahedron with an acute triangle 
facet opposite the apex vertex is a completely well-centered tetrahedron. 

Thus from any triangulation of a unit sphere S 2 with sufficiently small acute triangles we 
can create a completely well-centered tetrahedral mesh in M 3 by taking the cone u* r 2 of each 
acute triangle r 2 with the center of the sphere u. Figure 17 shows a completely well-centered 
tetrahedral mesh constructed in this fashion. The boundary of the mesh in Fig. 17 is an acute 
triangulation of S 2 selected from an infinite family of acute triangulations of S 2 . The next two 
paragraphs describe this family. 

Consider the set of vertices consisting of the north pole (0,0,1), the south pole (0,0,-1), 
and the vertices of two regular fc-gons, one in the plane z — 0.352 and the other in the plane 
z = —0.352. We set the polygons exactly off phase from each other. For instance, let the 
coordinates of the polygon vertices be 

/ 2iir\ / 2iir\ 

(0.936 cos — , 0.936 sin — ,0.352), i = 0, 1, . . . , k - 1 , 

(0.936 cos ( (2 '^ 1)7r ) ,0-936sin^ (2 ' + 1)7r j ,-0.352), i = 0, 1, . . . , k - 1 . 

Let each pole vertex be adjacent to all of the vertices of the closer regular polygon. This 
constructs fc isosceles triangles incident to each pole. We take each vertex of a regular polygon 
to be adjacent to the closer pole, the two neighbors on its own regular polygon, and two vertices 
from the other regular polygon. Triangles formed entirely from vertices of the two regular 
polygons are also isosceles. The example in Fig. 17 uses the result of this construction for the 
case k = 7. 

We claim that if fc > 4, then each triangle r 2 of this construction is acute and satisfies 
the condition that the distance from the origin to r 2 is greater than l/y/2. Since fc > 4 it 
is clear that the apex angles of the isosceles triangles incident to the poles are acute angles. 
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Verifying that the other triangles are acute and that the triangles are far enough from the origin 
is straightforward and we omit the details. Lemma 16 applies, and as an immediate consequence 
we have the following. 

Proposition 17. There are infinitely many triangulations of S 2 that can appear as the link of 
a vertex in a completely well-centered mesh. 

For large enough k, this construction of completely well-centered neighborhoods of a vertex 
using acute triangulations of a unit sphere S 2 can be generalized. One can use more than two 
regular fc-gons, alternating the phase between each successive fc-gon. 

We have seen that there are infinitely many triangulations of S 2 that cannot appear and 
infinitely many that can appear as the link of a vertex in a 3-well-centered mesh. The authors 
suspect that for m > 8 vertices the majority of triangulations of S 2 on m vertices are trian- 
gulations that can appear as a link of a vertex in a 3-well-centered mesh. We do not formally 
prove that conjecture in this paper, but in light of the the next proposition, it is highly likely; 
Proposition 18 provides a method for constructing new triangulations that can appear as the 
link of a vertex in a 3-well-centered tetrahedral mesh in W 3 . 

In Proposition 18 we consider a triangulation G of S 2 with a vertex of degree 3. In this 
context, the notation G — v\ refers to the triangulation of S 2 obtained by deleting vertex v\ and 
all faces incident to v\, replacing them with the face [w2 w 3 u 4], where v 2 ,v 3 ,V4 are the neighbors 
of v\ in G, ordered to keep the orientation consistent. 

Proposition 18. Let G be a triangulation of S 2 with a vertex V\ of degree 3, and let Vi, t>4 
be the neighbors of Vi in G. Let M be a tetrahedral mesh in R 3 consisting of a vertex u and its 
closed neighborhood Cl(St u), with Lk u isomorphic to G — V\. Lf 

(i) M is 3-well-centered 

(ii) face angle ZuViVj is acute for each i, j € {2,3,4}, i ^ j , 
then there exists a tetrahedral mesh M in K 3 and a vertex u of M such that 
(i) Lk u is isomorphic to G 
(ii) M is 3-well-centered 

(in) face angle ZuViVj is acute for each i,j <G {l,2,3,4},i ^ j . 

Proof. Figure 18 accompanies this proof and may help the reader understand the geometric 
constructions discussed in the proof. Consider a particular tetrahedral mesh that satisfies the 
conditions of the hypothesis. In this mesh the tetrahedron a = a 3 = [UV2V3V4] is 3-well-centered, 
so c(a) is interior to a. 

Let I be the line through u and c(a). Line £ intersects the circumsphere of a at two points. 
One of these is u, and the other we name u'. We define 

u' e = (1 — e)v! + eu , 

a point lying on I. Because a is 3-wcll-centercd, we know that segment uu' intersects triangle 
at some point u' £g , with 1/2 > e > 0. We can cut a into the three tetrahedra [uV2V3u' EQ ], 
[uv 3 V4u' So ], and [uVivW SQ \- 

For £0 > £ > we consider the three tetrahedra [uv2V3u' e ], [uv^v^u 1 ^, and [uv4V2u' e ]. We 
claim that for sufficiently small e > these three tetrahedra are 3-well-centered and the face 
angles Z.uu' e Vi, Luv{ii! e are acute for i — 2,3,4. 

Examining the face angles first, we note that at s = the circumcenters of the facets [ui^itg] 
coincide with c(a) and with each other. Indeed, each of these facets is a right triangle with 
its circumcenter lying on the hypotenuse uu' £ . As e increases, Zviuu' e does not change, ZuViU r £ 
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u' u' 

Figure 18: Given a 3- well-centered tetrahedron a = [UV2V3V4] with acute angles ZuViVj, one can 
construct three tetrahedra [UV2V3V1], [UV3V4V1], and [UV4V2V1] by adding a new vertex v\ = u' £ along the 
line £ through u and c(a). The circumcenters of the constructed tetrahedra lie along lines connecting 
c(cx) to the circumcenters c(rj) of the [uwj^j] facets of a. As discussed in Proposition 18, when v\ 
is close enough to v! — the reflection of u through c(a) — the constructed tetrahedra will be 3-well- 
centered and the angles Luv\Vi, ZuviVi will be acute. The angles ZviUVj do not need to be acute for 
this construction. For example, ZV2UV3 is not an acute angle in this figure. 



decreases, becoming smaller than ir/2, and Zuu' e Vi increases but remains less than tt/2 for 
sufficiently small e. 

Turning to the tetrahedra, then, we will argue that the specific tetrahedron [uv2V3u' £ \ is 
3- well-centered for sufficiently small e. An argument identical except for changed labels applies 
to the other two tetrahedra, so this will complete the proof. We know that, regardless of the 
value of e, the circumcenter of [uv2V3u' e ] lies on the line orthogonal to aff ([UV2V3]) = aff(T 4 ) 
passing through c(t4); this line is the locus of points equidistant from u, v 2 , and V3. The 
location of c([uv2V3u' £ ]) varies continuously with e. At e = 0, the circumcenter of tetrahedron 
[wu2U3Ug] coincides with c(<r), and as e increases from towards So, c{[uv2V3u' £ \) moves in the 
direction of vector c(ri) — c(a). Because ZUV2V3 and Zuw 3 w 2 are acute, we know that c(t 4 ) 
lies in the sector of aff(r4) interior to angle Zv 2 uv 3 . Thus segment c^cfa) n a is contained 
in [uv2V3u' £o ] C [uv2V 3 u' e ], and for sufficiently small e > 0, tetrahedron [uv2Vsu' e ] is 3-well- 
centered. □ 



Because the face angles A^uviVi, LuViV\ are acute in the construction of Proposition 18, the 
construction can be iterated. If a triangulation G of S 2 satisfies the conditions of Proposition 18, 
then a degree 3 vertex v\ can be inserted into face [W2W3W4]. In the new triangulation of S 2 , 
the three new faces incident to Vi satisfy the conditions of Proposition 18, so a degree 3 vertex 
can be inserted into any one of those three faces, and so on. In particular, starting from any 
completely well-centered mesh constructed from an acute triangulation of a unit sphere S 2 , one 
can successively insert vertices of degree 3 to create an infinite family of triangulations that can 
appear as the link of a vertex in a 3-well-centered mesh. 

It is also worth mentioning that each triangulation of a topological S 2 with 8 vertices 
V\, . . . ,vg that can appear as Lk u for a vertex u in a 3-well-centered mesh in M 3 has an 
embedding into M 3 for which all of the face angles AuviVj are acute for i, j G {1, . . . , 8}, i ^ j, 
where ViVj is an edge in the triangulation. Recall that there are 50 nonisomorphic triangulations 
of S 2 with 9 vertices [8]. Using Proposition 18 to add vertices of degree 3 to the various faces 
of triangulations of S 2 with 8 vertices, one can show that at least 34 of these 50 triangulations 
of S 2 with 9 vertices can appear as the link of a vertex in a 3-well-centered tetrahedral mesh 
embedded in R 3 . 
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6 Local Combinatorial Properties of 2-Well-Centered Tetra- 
hedral Meshes 



Corollary 14 shows that in a 3-well-centered mesh there are at least 7 edges incident to each 
vertex. In the following discussion we will see that the combinatorial constraints for a mesh to 
be 2-well-centered are quite different from the constraints for a mesh to be 3-well-centered, and 
in terms of the minimum number of edges incident to a vertex, are more stringent. As in Sec. 5, 
the discussion focuses on Lk u where u is a vertex interior to a tetrahedral mesh in R 3 . 

Definition. We say that a particular triangulation G of S 2 permits a 2-well-centered neighbor- 
hood of a vertex u if there exists a tetrahedral mesh M in M 3 such that u is an interior vertex 
of M, Lk u is isomorphic to G (as a simplicial complex) , and all facets of M incident to u are 
2-well-centered. (A facet means a 2-simplex in this context — a face of dimension n—1.) 

It should be noted that this definition does not directly address the question of whether the 
tetrahedra incident to u are 2-well-centered, since each tetrahedron incident to u has one facet 
lying on Lk u, and that facet is not incident to u. We shall see, however, that for tetrahedral 
meshes in M 3 , the smallest triangulation that permits a 2-well-centered neighborhood in the 
sense of this definition can, in fact, appear as the link of a vertex in a completely well-centered 
mesh. Finally, note that phrasing the problem in terms of the facets of M incident to u actually 
reduces the problem to determining whether the face angles at u are acute, because if there is 
an arrangement of rays at u such that all of the face angles formed at u by these rays are acute, 
then we can place the neighbors of u at the points where these rays intersect a unit sphere 
centered at u. This will create a neighborhood of u in which every 2-dimensional face incident 
to u is an isosceles triangle with an acute apex angle at u. 

The first result of this section is a simple observation that forms the foundation for the 
theory developed in the rest of the section. 

Lemma 19. Let u and v± be adjacent vertices in a tetrahedral mesh M embedded in K 3 and let 
Vi be a vertex of Lk u that is adjacent to v\. The angle /Lv\uvi is acute if and only if Vi G Hi, 
where Hi is the open halfspace that contains v\ and is bounded by the plane through u orthogonal 
to the vector v± — u. 

Proof. The angle Zviuvi is acute if and only if (vi — u,Vi — u) > 0, where (•, •) is the standard 
inner product on R 3 , and this holds if and only if Vi lies in Hi. □ 

The next two technical lemmas are based on Lemma 19. They lead to the proof of the 
main result of this section. In both lemmas and in the subsequent theorem we use the following 
notation. We denote by u a vertex in a tetrahedral mesh in M 3 , and the m vertices of Lk u 
are labeled vi, . . . , v m . For each vertex Vi, the plane through u orthogonal to Vi — u is denoted 
Pi, and the open halfspace bounded by Pi that contains Vi is denoted Hi. The other halfspace 
bounded by Pi will be called and we take this to be a closed halfspace, which contains its 
boundary Pi. The orthogonal projection of a vertex Vj into Pi will be denoted Pn(vj). 

Lemma 20. Let Vi and v 2 be nonadjacent vertices of Lk u, with V2 G H[. If Vi is a vertex 
adjacent to bothvi andv 2 such that Lviuvi and A.v 2 uVi are both acute angles, then the orthogonal 
projection of Vi into Pi lies in Pi n H 2 . 

Proof. The sketch in Fig. 19 illustrates this result. For an algebraic proof we assign a coordinate 
system with u as the origin and v\ lying on the positive z axis. Using coordinates (xi, yi, Zi) for 
vertex Vi, the condition v 2 € H[ means that z 2 < 0. Since the angle LviUVi is acute, Lemma 19 
implies that Vi must lie in Hi, and since the angle /Lv 2 uVi is acute, Lemma 19 implies that Vi 
must lie in H 2 . Thus Vi lies in Hi n H 2 . Since Hi n H 2 would be empty if v 2 had coordinates 
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Figure 19: If a 2-well-centered mesh contains two vertices v\ and v 2 that both lie in Lk u, are not 
adjacent to each other, and have a common neighbor v,- L and if v 2 lies in H[, then the orthogonal 
projection of Vi into Pi, i.e., the point Pr\(vi), must lie in Pi n H 2 . 




Figure 20: Let u be a vertex of a tetrahcdral mesh embedded in M 3 , and let vi, v 2 , vi, i>j be vertices 
of Lk u with adjacencies as shown. If the face angles at u between adjacent vertices of Lk u are all 
acute angles, but £viuv 2 is nonacute, then the projection of facet [viV 2 Vj] into Pi lies in Pi C\H 2 . 



(0,0, z 2 ), we can conclude that v 2 docs not lie on the z-axis. With the remaining freedom in 
defining a coordinate system we specify that v 2 has coordinates (x 2 ,0, z 2 ) with x 2 < 0. 

Now since Vi G Hi, we know that Zi > 0. We also know that (vi,v 2 ) = XiX 2 + ZiZ 2 > 0, 
because Vi G H 2 . We have established that ziz 2 < and that x 2 < 0. It follows that Xi < 0. The 
projection Pn(vi) has coordinates (a;,, y^, 0) and is interior to PiH H 2 — {(x, y,0) : x < 0}. □ 



Lemma 21. Let vi and v 2 be nonadjacent vertices of Lk u, with v 2 G H[. If [viV 2 Vj] is a 
2-simplex of Lk u, such that Vi,Vj are both adjacent to v\ and the face angles Lviuvi, ZviUVj, 
A.v 2 uVi, Z.v 2 uVj, are all acute angles, then Pri([viV 2 Vj]) C PiC\H 2 , i.e., the orthogonal projection 
of the entire facet [viV 2 Vj] into Pi lies in H 2 . 

Proof. See the sketch in Fig. 20. From the given hypotheses we can conclude by Lemma 20 that 
Pri(vi) and Pr\(vj) both lie in P\C\H 2 . Using the same coordinate system defined in the proof 
of Lemma 20, the point Pri(v 2 ) has coordinates (x 2 , 0, 0) with x 2 < 0, thus it lies in PiC\H 2 as 
well. It follows that the orthogonal projection of the facet [wiW2Vj] into Pi, which is the convex 
hull of Pri(vi), Pri(v 2 ), and Pri(vj), lies entirely in the convex set Pi(lH 2 . □ 



Applying the above two lemmas, we obtain a combinatorial necessary condition on the 
neighborhood of an interior vertex in a 2-well-centered mesh. 

Theorem 22. Let G be a triangulation of S 2 with m vertices. If G contains a vertex v\ of 
degree d(vi) > m — 3, then G does not permit a 2-well-centered neighborhood. 
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Figure 21: When Lk u has m vertices and one of the vertices «i has degree d(v\) > m — 3, any 
geometric realization of Cl(St u) in M 3 with all face angles at u acute is not an embedding. Theorem 22 
shows that if we consider such a geometric realization and project every facet that intersects H[ into 
Pi, then the union of the projected facets does not contain u. The sketch at left shows an example of 
a geometric realization of a tetrahedral mesh Cl(St u) in M 3 such that every face angle at vertex u is 
acute. In the sketch, v\ has degree d(v\) = 6 = to — 3 in Lk u. The sketch at right shows the result 
of taking the geometric realization on the left and projecting each facet that intersects H[ into P\. 

Proof. We consider a vertex u such that Lk u is isomorphic to G where G has a vertex of degree 
at least to — 3 and consider a geometric realization of Cl(St u) in M? . Label the vertices of Lk u 
with the labels v±, v 2 , • • • , v m such that V\ is a vertex of maximum degree and the (at most two) 
vertices not adjacent to v\ are listed immediately after v\ (e.g., labeled v 2 , «3 if there are two 
of them). We choose a coordinate system on K 3 such that u is at the origin and vi lies on the 
positive z axis. 

Assume that all of the face angles ZviUVj are acute. We claim this implies that for any 
facet [viVjVk] with at least one vertex in H[, the orthogonal projection of the facet into Pi, 
i.e., Pri([viVjVk]), does not contain vertex u. Assuming this claim for the moment, we see that 
u lies outside the (solid) polyhedron bounded by Lk u. (See Fig. 21.) Since u is outside this 
polyhedron, some 3-simplex incident to u must be inverted. Thus the geometric realization of 
Cl(St u) is not an embedding, and the claim completes the proof. 

We proceed to prove the claim. Noting that v\ € Hi by our definition of Hi , we observe that 
for i > 4, vertex Vi must lie in Hi because Vi is adjacent to v\. (This follows from Lemma 19.) 
Thus there are only two types of facets that may have nonempty intersection with H[. The 
first type is [v^Vj] or [i>ji>3Vj] where Vi and Vj both are adjacent to vi, and the second type is 
[w 2 W3«j] for j > 4. Consider, then, the first type of facet, taking the specific notation [viV 2 Vj]. 
(The same argument applies to [viV^Vj].) If v 2 lies in Hi, we are done; the facet does not 
intersect H[. Otherwise v 2 lies in H[. Hence Zviuv 2 is nonacute, and v 2 is not adjacent to v\. 
Lemma 21 applies. 

The proof for facets of the second type is more complicated. If both v 2 and t>3 lie in H±, 
we are done. If one vertex lies in Hi U Pi and the other lies in H[, we assume without loss of 
generality that z 2 < and z 3 > 0. 

Then v 2 is not adjacent to vi. If V3 is adjacent to v\, then Lemma 21 applies directly with 
t>3 functioning as t^. On the other hand, even if v% is not adjacent to v\, the arguments of 
Lemmas 20 and 21 can be applied with w 3 functioning as Vi. (In the proofs of Lemmas20 and 21 
we used Vi adjacent to vi to establish only that > and that v 2 does not lie on the z-axis. 
The latter holds in this case because v 2 and vi have common neighbor Vj ^ V3.) 

This leaves the case z 2 < and z 3 < 0. As noted above, v 2 does not lie on the z-axis. We 
choose the coordinate system with v 2 = (22,0,2:2); x 2 < 0. We also assume without loss of 
generality that j/3 > 0. (We can reflect through the plane y = if 2/3 < 0.) See Fig. 22 for 
sketches related to this case. 
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Figure 22: In the proof of Theorem 22, the most difficult case to analyze is a facet of the second type 
[i>2t>3fj] with Z2 < and z 3 < 0. These sketches illustrate the projection of the facet [i>2i>3Uj] m t° -Pi 
for the two subcases x 3 < (left) and x 3 > (right). In both subcases, Pri(vj) lies in Pi n H 2 n -ff 3 . 
When X3 < 0, Pri(v 3 ) E Pi H2, so the projection of facet [f2^3«j] into Pi is a subset of Pi (~l i?2- 
When X3 > 0, the projected facet can be decomposed into two pieces meeting along \poPri(vj)]. One 
piece lies in Pi(lH 2 , and the other piece lies in Pi(lH 3 . In both subcases we see that u £ Pri([v2V 3 Vj]). 

By applying Lemma 19 three times, we obtain Vj € Hif]H2^H 3 , and by applying Lemma 20 
twice we obtain Pr\(vj) e Pi n H 2 n H 3 . If x 3 < 0, then the whole segment Pri([v 2 v 3 ]) lies in 
PiC\H 2 , and since Pn(vj) e PiC\H 2 , it follows that Pri([u 2 «3«j]) C PiC\H 2 . Thus Pri([t) 2 «3Uj]) 
does not contain u. 

So we assume that a: 3 > 0. Now if x 3 = we know that 1/3^0 because v 3 , like w 2 , does 
not lie on the z-axis. Moreover, x 3 > also implies y 3 ^ 0, since otherwise Pi n 7?3 would be 
{(x, y, 0) : x > 0}, yielding P 1 n i? 2 n i? 3 = 9 Pri(t)j). A point on Pri([w 2 w 3 ]) has the form 



with < A < 1. Thus for a point p = (x p , y p , 0) on Pri([v2V 3 ]), either x p < and the point lies 
in Pi n H 2 or both x p > and y p > so that (p, w 3 ) = x p a;3 + y p y 3 > and the point lies in 
Pi n H 3 . We conclude that Pri([v 2 v 3 ]) C Pi n (H 2 U H 3 ). 

Finally we note that there must exist a point po — (e, y Po ,0) with e < such that po lies on 
Pri([t)2^3]) and po G Pi fl P 2 H P 3 . Thus we can decompose Pri([v2V 3 Vj]) into the two pieces 
Po * P r i([ v 2Vj]) and p * Pri([w3^']), with the first piece lying in Pi n P2 and the second piece 
lying in Pi n H 3 . It follows that Pri([v2V 3 Vj]) does not contain u. □ 

Recall that Euler's formula specifies a relationship between the number of vertices, edges, and 
faces in a planar graph. If m, e, and / are the number of vertices, edges, and faces respectively, 
then Euler's formula states that m — e + f = 2 for planar graphs. In a planar triangulation 
each face is incident to three edges and each edge is incident to two faces, so 2e = 3/, and 
the relationship / = 2(m — 2) can be derived. Moreover, in a planar triangulation each face is 
incident to three vertices, and vertex Vi is incident to d(vi) faces, so ^ d(vi) = 3/ = 6(m — 2). 
Combining these consequences of Euler's formula with Theorem 22, we easily obtain a lower 
bound on the number of edges incident to an interior vertex in a 2-well-centered tetrahedral 
mesh in R 3 . 

Corollary 23. Let M be a 2-well-centered (or completely well- centered) tetrahedral mesh in M 3 . 
For every vertex u interior to M, at least 9 edges of M are incident to u. 

Proof. Let G — Lk u for some interior vertex u of a 2-well-centered mesh M, and let m be the 
number of edges incident to u, i.e., the number of vertices of G. Consider the possibility m = 8. 



XPri(v 2 ) + (1 - X)Pr 1 (v 3 ) = A(x 2 , 0, 0) + (1 - A) (x 3 , y 3 , 0) 
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Figure 23: A completely wcll-ccntcrcd mesh with an interior vertex u such that Lk u has nine 
vertices and degree list (5,5,5,5,5,5,4,4,4). The vertex coordinates are listed in the table at 
right; vertex u is at the origin. 



Euler's formula shows that for ra = 8we have J2i d(vi) = 36 so the average vertex degree is 4.5, 
and there must be at least one vertex of degree at least 5 = m — 3. By Theorem 22, this cannot 
occur, for such a graph G would not permit a 2- well-centered neighborhood of u. Similarly, if 
m = 7 the average degree is 30/7 > 4 and there must be a vertex of degree at least m — 2. For 
m = 6 the average degree is 4 and there must be a vertex of degree at least m — 2. In each of 
the cases m = 5 and m = 4, there is only one triangulation, and this triangulation has a vertex 
of degree m — 1 . □ 



When m = 9, the average degree is 4|, and there is a triangulation of S 2 with degree list 
(5,5,5,5,5,5,4,4,4) that permits a completely well-centered neighborhood. Figure 23 shows a 
figure of a completely-well-centered mesh that has a single interior vertex u such that Lk u is a 
9- vertex triangulation of S 2 with the specified degree list. 

We have already seen that there are infinitely many triangulations of S 2 that can appear as 
the link of an interior vertex in a 2- well-centered mesh (Proposition 17). In the spirit of Proposi- 
tion 18, we now discuss some ways to use an existing triangulation that permits a 2- well-centered 
neighborhood to construct new triangulations that permit a 2-wcll-centcrcd neighborhood. The 
next two propositions show that one can add vertices of degree 3, subtract vertices of degree 3, 
or add vertices of degree 4 to obtain new triangulations that permit a 2-well-centered neighbor- 
hood. In Proposition 24 we again use the notation G — v\ used in Proposition 18. 

Proposition 24. A triangulation G of S 2 that contains a vertex v\ of degree three permits a 
2-well-centered neighborhood if and only if the triangulation G — v\ permits a 2-well-centered 
neighborhood. 

Proof. First we suppose that G permits a 2-well-centcrcd neighborhood. Then consider some 
tetrahedral mesh embedded in IR 3 that contains a vertex u with Lk u isomorphic to G and all 
face angles ZviUVj acute. We choose a coordinate system on R 3 such that u lies at the origin 
and identify each vertex Vi of Lk u with the vector originating at the origin and terminating at 
Vi. Now vector v\ makes an acute face angle for each of the three facets that are incident to the 
edge [uvi). Deleting v\ from Lk u removes the three facets that are incident to edge [uvi], but 
has no effect on the other facets incident to u or face angles at u. Thus all facets incident to u 
remain acute after removing v\ , and the modified neighborhood of u is a mesh embedded in R 3 
that certifies that G — vi permits a 2-well-centered neighborhood. 
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On the other hand, if we suppose that G — v\ permits a 2-wcll-centcrcd neighborhood, we 
will be able to add vertex v\ and still have all face angles at u acute. Consider some specific 
tetrahedral mesh embedded in R 3 containing a vertex u such that Lk u is isomorphic to G — v\ . 
Let v 2 ,v 3 , and v 4 be the three vertices of Lk u that are adjacent to v\ in G. Then the mesh 
contains facets [W2UU3], [v 3 uv 4 ], and [w 4 uw 2 ]- Moreover, since the face angles at u are acute, 
we have (vi,Vj) > for each € {2,3,4} x {2,3,4}. It follows that if we insert vertex v\ 

satisfying v\ = X 2 v 2 + X 3 v 3 + X4V4 with each A, > 0, then for i = 2, 3, 4 we have 

(vi,Vi) = X 2 (v2,Vi) + X 3 (v 3 ,Vi) + A 4 (u 4 ,w i ) > 0. 

In other words, as long as v\ lies interior to the cone at u bounded by vectors v 2 , t>3, and w 4 , it 
will make acute face angles with each of v 2 , v 3 , and V4. □ 



Notice that Proposition 24 also implies that adding or deleting a degree three vertex from 
a triangulation of S 2 that does not permit a 2-well-centered neighborhood creates another 
triangulation of S 2 that does not permit a 2-well-centered neighborhood. In particular, this 
means that Theorem 22 does not characterize the triangulations of S 2 that cannot appear as 
the link of a vertex in a 2-well-centered tetrahedral mesh in M 3 . 

In the next proposition, we consider the case of a triangulation G of S 2 with a vertex v\ 
such that d(v\) = 4. To talk about removing vertex v\ from G in this case, we need to specify 
an edge to add after removing the vertex. Let v 2 , v 3 , w 4 , and V5 be the neighbors of v\, listing 
in cyclic order. Then (G — V\) U [^2^4] is the triangulation of S 2 obtained from G by removing 
vertex v\ along with the four edges and triangles incident to v\ and adding the edge [V2V4] along 
with the two triangles [V2U3V4] and [t^^^s]- 

Proposition 25. Consider a triangulation G of S 2 that contains a vertex V\ of degree four 
with neighbors v 2 , V3, V4, v$ (listed in clockwise order). If (G — Vi) U [«2« 4 ] or (G — U [V3W5] 
permits a 2-well-centered neighborhood, then G permits a 2-well-centered neighborhood. 

Proof. Suppose without loss of generality that (G — vi) U [v 2 V4] permits a 2-well-centered neigh- 
borhood. Let u be a vertex for which Lk u is isomorphic to (G — v\) U [^2^4] and consider 
some embedding of u * Lk u into M 3 such that all face angles at u are acute. We choose a 
coordinate system on M 3 such that u lies at the origin and identify each vertex of Lk u with the 
vector originating at the origin u and terminating at the vertex. We know that (^2,^3) > 0, 
(«3,v 4 ) > 0, (u4,«5) > 0, (v5,v 2 ) > 0, and (1^2,^4) > 0, because each pair of vectors bounds a 
face with an acute face angle at u. 

Now let v\ = (v 2 + V4)/2 and, deleting the facet [v2«v 4 ], add the four facets for 
i = 2, 3, 4, 5. The new facets [vi W2] and [viwv 4 ] have face angles at u that are smaller than the 
face angle Lv 2 uv4 was, so they are acute. The facets [U1WU3] and [uiuws] also have acute face 
angles at u because 

(vi,v 3 ) = ~(v 2 ,v 3 ) + -(v4,v 3 ) > 

and similarly 

( v i, v 5) = ^(v 2 ,v 5 ) + ~(v4,v 3 ) > 0. 

We see that adding v\ = (v 2 + w 4 ) /2 has created a new mesh for which all face angles at u are 
acute. Thus G permits a 2-well-centered neighborhood. □ 
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Figure 24: No 3-well-centered mesh of 
the cube has a face with this triangula- 
tion. 




Figure 25: No 2-well-centered or 3-well- 
centercd mesh of the solid cube has a 
tetrahedron with three of its facets lying 
in faces of the cube. 



7 Applications to the Cube 

The theoretical results presented in this paper are useful for creating well-centered meshes of 
specific regions in W 3 . In particular, one might design a tetrahedral mesh of a volume so that 
it meets all of the combinatorial conditions discussed in Sees. 5 and 6. Then applying the 
optimization procedure discussed in [17], one may hope to obtain a well-centered mesh of the 
domain. This technique was successfully used in [15] to create well-centered meshes of several 
domains in M 3 , including the cube. 

The theory developed in this paper has several obvious implications for the combinatorial 
properties of a well-centered triangulation of the cube. For example, no cube corner tetrahedron, 
e.g., the tetrahedron shown in Fig. 25, can be 3-well-centered; considering the bottom facet to 
be a given facet, we see that the fourth vertex of the tetrahedron projects onto (not inside) the 
circumcircle of the given facet, violating the necessary Cylinder Condition of Proposition 3. 

It follows that in a 3-well-centered mesh of the cube there must be at least two tetrahedra 
incident to each corner of the cube. Indeed, there must be at least three tetrahedra incident to 
each corner of the cube in a 3-well-centered mesh. In the case of two tetrahedra incident to a 
corner vertex there must be exactly four edges incident to the corner vertex, of which three are in 
the directions of the coordinate axes. The fourth edge must lie in a face, and both tetrahedra are 
incident to the axis orthogonal to the face containing the fourth edge. The Cylinder Condition 
applies again, and we see that the mesh cannot be 3-well-centered. 

Ad hoc arguments from basic Euclidean geometry provide more restrictions on well-centered 
triangulations of the solid cube. For instance, in any 3-well-centered mesh of the cube, no face 
of the cube is triangulated as shown in Fig. 24, with two right triangles meeting along the 
hypotenuse. The two right triangles have the same circumccntcr, which lies at the midpoint 
of the common hypotenuse of the triangles the center of the face of the cube. For cither 
triangle, a tetrahedron having that triangle as a facet must have its circumcenter on a line I 
perpendicular to the face of the cube that meets the cube face at its center. Considering two 
tetrahedra o\ and a 2 , each having one of the right triangles as a face, it can be shown that at 
most one of U\ , a 2 can be 3-well-centered. There is a plane that contains the hypotenuse of the 
right triangles and divides M 3 into two open halfspaces Hi and H 2 such that ci G Cl(ifi) and 
u 2 G Cl(H 2 ). If ci and a 2 share a common face, the plane must be aff(cri n 02), but otherwise 
there is some flexibility in the choice of the plane. The portion of I interior to the cube is 
cither in the boundary between Hi and H 2 or without loss of generality can be assumed to lie 
entirely in Hi. In cither case, the circumcenter of cr 2 is not strictly interior to a 2 , so a 2 is not 
3-well-centered. 

The paper [15], along with discussing some well-centered meshes of the cube, raises the 
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question of the minimum number of tetrahedra needed to create a well-centered triangulation 
of the cube. We can use the statements above to derive some simple lower bounds on the 
number of tetrahedra in a well-centered mesh of the 3-cube. In a triangulation of the cube, the 
number of tetrahedra incident to the surface of the cube is a lower bound on the total number 
of tetrahedra, so one can obtain a lower bound on the number of tetrahedra by counting the 
number of triangular facets in a surface triangulation. The number of facets is not a direct 
lower bound, since there may be a single tetrahedron with multiple facets in the surface of 
the cube. Because there are at least three distinct tetrahedra incident to a cube corner in a 
3-well-centered triangulation of the cube, a tetrahedron cannot be counted more than twice in 
counting the number of surface facets of a 3-well-centered triangulation of the cube. The same 
holds true for 2-well-centered triangulations of the cube, since three of the facets of a cube corner 
tetrahedron are right triangles. 

Noting, then, that each face of the cube must contain at least 3 triangles in a 3-well-centered 
mesh of the cube and at least 8 triangles in a 2-well-centered mesh of the cube, we easily obtain a 
lower bound of 9 tetrahedra for a 3-well-centered triangulation of the cube, and 24 tetrahedra for 
a 2-well-centered triangulation of the cube. (These lower bounds arc mentioned in [15] without 
the details of the geometric or combinatorial arguments.) It should be possible to improve both 
of these bounds, but these relatively simple bounds help demonstrate a possible application of 
this paper's theory and are a starting place for a more careful analysis. 

8 Conclusions 

In this paper we introduced several geometric propositions related to n-well-centered simplices 
and gave an algebraic characterization of an n-well-centered simplex in terms of cubic polynomial 
inequalities. We applied the geometric propositions to the study of the combinatorial properties 
of well-centered meshes, especially well-centered tetrahedral meshes. 

We considered triangulations of topological S 2 and showed that the set of such triangulations 
that cannot appear as the link of a vertex in a 3-wcll-ccntcrcd (or 2-wcll-centered or completely 
well-centered) tetrahedral mesh embedded in K 3 is an infinite set, contrasting this to the anal- 
ogous question for triangle meshes in R 2 . We showed also that the set of triangulations of S 2 
that do appear as the link of a vertex in some completely well-centered (or 3-well-centered or 
2-well-centered) tetrahedral mesh embedded in R 3 is an infinite set. We proved several results 
in the direction of classifying which triangulations of S 2 can appear as the link of a vertex in a 
2-well-centered or 3-well-centered tetrahedral mesh embedded in R 3 . 

The work on combinatorial properties of well-centered meshes leads to some interesting open 
questions. Is there a compact way to express a complete characterization of which triangulations 
of S 2 can appear in a 3-well-centered (or 2-wcll-centered or completely well-centered) mesh 
embedded in IR 3 ? Is the necessary condition described in Theorem 13 a complete characterization 
for vertex links in 3- well-centered tetrahedral meshes in R 3 ? If a triangulation of S 2 permits 
a 2-well-centered neighborhood in the sense defined in this paper, does this imply that it can 
appear as the link of a vertex in a 2-well-centered tetrahedral mesh embedded in K 3 ? If a 
triangulation of S 2 can appear as the link of a vertex in both a 2-well-centered tetrahedral mesh 
in R 3 and a 3-well-centered tetrahedral mesh in R 3 , does this guarantee that it can appear as 
the link of a vertex in a completely well-centered mesh? 

Beyond the questions about tetrahedral meshes there are questions about higher dimensions. 
Is it possible to extend the results of Sec. 6 to say something about 2-well-centered meshes in 
higher dimensions? Certainly Lemmas 19, 20, and 21 can be generalized to higher dimensions. Is 
it the case that for each n there exists a completely well-centered n-simplicial neighborhood of a 
vertex embedded in R n ? If so, Lemma 16 may provide a way of constructing such neighborhoods. 
We note, however, that for n > 5 there is no dihedral acute n-simplicial neighborhood of a vertex 
embedded in R n [6] [7]. 

This work also leads to some interesting practical questions about creating well-centered 
tetrahedral meshes in R 3 . In particular, now that there is some understanding of which tri- 
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angulations can appear as the link of a vertex in a well-centered mesh, we hope that practical 
methods can be developed for improving the local mesh connectivity of tetrahcdral meshes. An 
algorithm for improving mesh connectivity of triangle meshes in R 2 appears in [16], but it is 
not obvious how to formulate that type of algorithm for tetrahcdral meshes in R 3 . It is also 
worth noting that a triangulation of S 2 that can theoretically appear as the link of a vertex 
in a well-centered mesh might be a poor neighborhood for a vertex in a practical setting. For 
instance, in a triangle mesh in R 2 , a vertex link with 100 vertices can appear in a 2- well-centered 
triangle mesh embedded in R 2 , but a mesh containing such a vertex link would be considered 
poor quality in almost any application. Is there a good way to rate vertex links in tetrahedral 
meshes according to their applicability in a practical setting? 

We hope that our results will motivate others to investigate these interesting and important 
questions. 
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